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l/"") ' Abstract 

The purpose of this paper is to develop a synthesis theory for linear dynamical quantum 
stochastic systems that are encountered in linear quantum optics and in phenomcnological 
models of linear quantum circuits. In particular, such a theory will enable the systematic 
realization of coherent /fully quantum linear stochastic controllers for quantum control, 
amongst other potential applications. We show how general linear dynamical quantum 
stochastic systems can be constructed by assembling an appropriate interconnection of one 
^ ' degree of freedom open quantum harmonic oscillators and, in the quantum optics setting, 

O" 1 discuss how such a network of oscillators can be approximately synthesized or implemented 

in a systematic way from some linear and nonlinear quantum optical elements. An example 
£Nj ■ is also provided to illustrate the theory. 
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O ■ 1 Background and motivation 



In recent years there has been an explosion of interest in exploitation of quantum mechanical 
systems as a basis for new quantum technologies, giving birth to the field of quantum infor- 
mation science. To develop quantum technologies, it has been recognized from early on that 
quantum control systems will play a crucial role for tasks such as manipulating a quantum me- 
chanical system to perform a desired function or to protect it from external disturbances [TJ [2] . 
Moreover, recent advances in quantum and nanotechnology have provided a great impetus for 
research in the area of quantum feedback control systems; e.g., see O HI [5j [H [8] . 

Perhaps just about the simplest and most tractable controller to design would be the lin- 
ear quantum controllers, and this makes them an especially attractive class of controllers. In 
this class, one can have classical linear quantum controllers that process only classical signals 
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which are obtained from a quantum plant by measurement of some plant output signals (e.g., 
[5j[9l[10]), but more recently there has also been interest in fully quantum and mixed quantum- 
classical linear controllers that are able to manipulate quantum signals |11] fT2| [T3| W$\ [T5] . In 
fact, an experimental realization of a fully quantum controller in quantum optics has been 
successfully demonstrated in [8]. As noted in that paper, the class of fully quantum controllers 
or coherent-feedback controllers, as they are often known in the physics literature, presents 
genuinely new control-theoretic challenges for quantum controller design. An important open 
problem raised in the works [12|, [T3l [Til [T5] is how one would systematically build or imple- 
ment a general, arbitrarily complex, linear quantum controller, at least approximately, from 
basic quantum devices, such as quantum optical devices. This problem can be viewed as a 
quantum analogue of the synthesis problem of classical electrical networks (in this paper the 
qualifier "classical" refers broadly to systems that are not quantum mechanical) that asks the 
question of how to build arbitrarily complex linear electrical circuits from elementary passive 
and active electrical components such as resistors, capacitors, inductors, transistors, op-amps, 
etc. Therefore, the quantum synthesis problem is not only of interest for the construction of 
linear quantum stochastic controllers, but also more broadly as a fundamental aspect of linear 
quantum circuit theory that arises, for example, in quantum optics and when working with 
phenomenological models of quantum RLC circuits such as described in [16], as well as in 
relatively new fields such as nanomechanical circuit quantum electrodynamics \17\ [15] . 

A key result of this paper is a new synthesis theorem (Theorem[2|) that prescribes how an ar- 
bitrarily complex linear quantum stochastic system can be decomposed into an interconnection 
of basic building blocks of one degree of freedom open quantum harmonic oscillators and thus 
be systematically constructed from these building blocks. In the context of quantum optics, we 
then propose physical schemes for "wiring up" one degree of freedom open quantum harmonic 
oscillators and the interconnections between them that are required to build a desired linear 
quantum stochastic system, using basic quantum optical components such as optical cavities, 
beam splitters, squeezers, etc. An explicit yet simple example that illustrates the application 
of the theorem to the synthesis of a two degrees of freedom open quantum harmonic oscillator 
is provided. 

1.1 Elements of linear electrical network synthesis 

To motivate synthesis theory in the context of linear dynamical quantum systems, we start 
with a brief overview of aspects of linear electrical network synthesis that are relevant for the 
current work. 

As is well known, a classical (continuous time, causal, linear time invariant) electrical net- 
work described by a set of (coupled) ordinary differential equations can be analyzed using vari- 
ous representations, for example, with a frequency domain or transfer function representation, 
with a modern state space representation and, more recently, with a behavioral representation. 
It is well known that the transfer function and state space representation are equivalent in 
the sense that one can switch between one representation to the other for any given network. 
However, although one can associate a unique transfer function representation to a state space 
representation, the converse is not true: for a given transfer function there are infinitely many 
state space representations. The state space representation can be made to be unique (up to 
a similarity transformation of the state space matrices) by requiring that the representation 
be of minimal order (i.e., the representation is both controllable and observable). The syn- 
thesis question in linear electrical networks theory deals with the inverse scenario, where one 
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Figure 1: Schematic for the implementation of the classical system ([1} 




Figure 2: Hardware implementation of the schematic diagram shown in Figured) 



is presented with a transfer function or state space description of a linear system and would 
like to synthesize or build such a system from various linear electrical components such as 
resistors, capacitors, inductors, op-amps, etc. A particularly advantageous feature of the state 
space representation, since it is given by a set of first order ordinary differential equations, is 
that it can be inferred directly from the representation how the system can be systematically 
synthesized. For example, consider the system below, given in a state space representation: 



dx(t) 
dt 

y(t) 



2 5 
-2 -4 



x(t) + 
[01] x(t)+u(t), 



1 

0.1 



u(t), 



(1) 



where x(t) is the state, u(t) is the input signal, and y(t) is the output signal. In an electrical 
circuit, u{t) could be the voltage at certain input ports of the circuit and y(t) could be the 
voltage at another set of ports of the circuit, different from the input ports. This system can be 
implemented according to the schematic shown in Figure [H This schematic can then be used to 
to implement the system at the hardware level as shown in Figure [21 [19\ Chapter 13]. However, 
linear electrical network synthesis is a mature subject that deals with much more than just 
how one can obtain some realization of a particular system. For instance, it also addresses 
fundamental issues such as how a passive network, a network that does not require an external 
source of energy, can also be synthesized using only passive electrical components, and how to 
synthesize a given circuit with a minimal number of circuit elements or with a minimal number 
of certain types of elements (such as active elements). In this paper our primary objective is to 
develop an analogously systematic method for synthesizing arbitrarily complex linear dynamical 
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quantum stochastic systems that are given in an abstract description that is similar in form 
to dl]). These linear dynamical quantum stochastic systems are ubiquitous in linear quantum 
optics, where they arise as idealized models for linear open quantum systems. However, since 
there is currently no comprehensive synthesis theory available for linear dynamical quantum 
systems (as opposed to static linear quantum systems in linear quantum optics that have been 
studied in, e.g., [20J) and related notions such as passivity have not been extensively studied 
and developed, here we focus our attention solely on the development of a general synthesis 
method that applies to arbitrary linear dynamical quantum systems which does not exploit 
specific physical properties or characteristics that a particular system may possess (say, for 
instance, passivity). Although the latter will be an important issue to be dealt with in further 
development of the general theory, it is beyond the scope of the present paper (which simply 
demonstrates the existence of some physical realization). 

1.2 Open quantum systems and quantum Markov models 

A quantum system is never completely isolated from its environment and can thus interact 
with it. Such quantum systems are said to be open quantum systems and are important in 
modeling various important physical phenomena such as the decay of the energy of an atom. 
The environment is modeled as a separate quantum system in itself and can be viewed as a heat 
bath to which an open quantum system can dissipate energy or from which it can gain energy 
(see [21\ Chapters 3 and 7] ) . An idealization often employed in modeling the interaction between 
an open quantum system and an external heat bath is the introduction of a Markovian assump- 
tion: the dynamics of the coupled system and bath is essentially "memoryless" in the sense that 
future evolution of the dynamics of the coupled system depends only on its present state and 
not at all on its past states. Open quantum systems with such a property are said to be Markov. 
The Markov assumption is approximately valid under some physical assumptions made on the 
system and bath, such as that the heat bath is so much "larger" than the system (in the sense 
that it has many more degrees of freedom than the system) and is weakly coupled to the system 
that its interaction with the latter has little effect on its own dynamics and can thus be ne- 
glected; for details on the physical basis for this Markovian assumption, see |2 1 1. Chapters 3 and 
5]. Markov open quantum systems are important, as they are often employed as very good ap- 
proximations to various practically relevant open quantum systems, particularly those that are 
encountered in the field of quantum optics, yet at the same time are relatively more tractable to 
analyze as their dynamics can be written in terms of first order operator differential equations. 

In Markov open quantum systems, heat baths can be idealistically modeled as a collection 
of a continuum of harmonic oscillators oscillating at frequencies in a continuum of values. An 
important consequence of the Markov approximation in this model is that the heat bath can be 
effectively treated in a quantum statistical sense as quantum noise [211 section 3.3], and thus 
Markov open quantum systems have inherently stochastic quantum dynamics that are most 
appropriately described by quantum stochastic differential equations (QSDE) [2T | [22 | [23" 1 [24"]. 
To be concrete, a single heat bath in the Markov approximation is formally modeled as an 
operator- valued quantum white noise process n(t), where t > denotes time, that satisfies 
the singular commutation relation [n(t),r](t)*] = 5(t — t'), where * denotes the adjoint of an 
operator, 8(t) is the Dirac delta function, and the commutator bracket [•, •] acts on operators A 
and B as [A, B] = AB — BA. Examples of heat baths that have been effectively modeled in such 
a way include vacuum noise, squeezed and laser fields in quantum optics [21], and infinitely long 
bosonic transmission lines [16] . See also [25] for a brief intuitive overview of the modeling of a 
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free-traveling quantized electromagnetic wave as quantum white noise. The formal treatment 
with quantum white noises can be made mathematically rigorous by considering the bosonic 
annihilation process A(t) (on a Fock space) that can be formally defined as the "integral" of 
rj(t): A(t) = J"q7/(s)cZs and its adjoint process A*(t) = A(t)*. We shall refer to the operator 
process A(t) simply as a bosonic field. The celebrated Hudson-Parthasarathy (H-P) stochastic 
calculus provides a framework for working with differential equations involving the processes 
Ait) and A*(t), as well as another fundamental process on a Fock space called the gauge 
process, denoted by A(t), that models scattering of the photons of the bosonic heat bath (at a 
formal level, one could write A(t) = f r](s)*r](s)dt). More generally, a quantum system can be 
coupled to several independent bosonic fields A\(t), . . . , A n (t), with Aj(t) = f rjj(s)ds, and in 
this case there can be scattering between different fields modeled by interfield gauge processes 
Ajk(t) = Jq Vj( s )*Vk{s)ds (in interfield scattering, a photon is annihilated in one field and then 
created in another). 

1.3 Linear dynamical quantum stochastic systems 

Linear dynamical quantum stochastic systems (e.g., see [101 112] ) arise in practice as idealized 
models of open quantum harmonic oscillators whose canonical position and momentum oper- 
ators are linearly coupled to one or more external (quantum) heat baths (the mathematical 
modeling involved is discussed in section [2]). Here a quantum harmonic oscillator is a quantized 
version of a classical harmonic oscillator in which the classical position and momentum vari- 
ables q c and p c , respectively, are replaced by operators q, p on an appropriate Hilbert space (in 
this case the space L 2 (1R)) satisfying the canonical commutation relations (CCR) [q,p] = 2i. It 
is said to be open if it is interacting with elements of its environment. For instance, consider 
the scenario in [5] of an atom trapped in an optical cavity. The light in the cavity is strongly 
coupled to the atomic dipole, and as the atom absorbs and emits light, there are random me- 
chanical forces on the atom. In an appropriate parameter regime, the details of the optical and 
atomic dipole dynamics are unimportant, and the optical field can be modeled as an environ- 
ment for the atomic motion. Under the assumptions of [5] the "motional observables" of the 
trapped atom (its position and momentum operators) can then be treated like those of an open 
quantum harmonic oscillator. Linear Markov open quantum models are extensively employed 
in various branches of physics in which the Markov type of arguments and approximations 
such as discussed in the preceding subsection can be justified. They are particularly prominent 
in quantum optics, but have also been used, among others, in phenomenological modeling of 
quantum RLC circuits [16], in which the dissipative heat baths are realized by infinitely long 
transmission lines attached to a circuit. For this reason, the general synthesis results developed 
herein (cf. Theorem [2|) are anticipated to be be relevant in various branches of quantum physics 
that employ linear Markov models. For example, it has the potential of playing an important 
role in the systematic and practical design of complex linear photonic circuits as the technology 
becomes feasible. 

A general linear dynamical quantum stochastic system is simply a many degrees of freedom 
open quantum harmonic oscillator with several pairs of canonical position and momentum 
operators (?fc,Pfc> with k ranging from 1 to n, where n is the number of degrees of freedom of the 
system, satisfying the (many degrees of freedom) CCR [qj,Pk] = 2io~jk and [qj,qk] = \Pj-,Pk\ = 0, 
where 8^ is the Kronecker delta which takes on the value unless j = k, in which case it takes 
on the value 1, that is linearly coupled to a number of external bosonic fields Ai, . . . ,A m . In 
the interaction picture with respect to the field and oscillator dynamics, the operators qj,Pj 
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A(t) 



*y(t) 

Figure 3: A Fabry-Perot optical cavity. The black rectangle denotes a mirror fully reflecting at 
the cavity resonance frequency, while the white rectangle denotes a mirror partially transmitting 
at that frequency. 




evolve unitarily in time as qj(t),pj(t) while preserving the CCR [qj(t),Pk(t)] = 2i°~jk an d 
Qk(t)] = \Pj{t)iPk(t)] — Vt > 0, and the dynamics of the oscillator is given by (here 
x(t) = ( qi (t), Pl (t),...,q n (t),p n (t)) T and A(t) = (A^t), . . . , A m (t)) T ) 



dx(t) = Ax(t)dt + B 



dA(t) 
dA(t)* 

dy(t) = Cx{t)dt + DdA(t), (2) 



where A G R nxn , B G C nx2m , C G C mxn , and D G C mxm . Here the variable y(t) acts as the 
output of the system due to interaction of the bosonic fields with the oscillator; a component 
Uj(t) of y(t) is the transformed version of the field Aj(t) that results after it interacts with 
the oscillator. Hence, Aj(t) can be viewed as an incoming or input field, while yj(t) is the 
corresponding outgoing or output field. To make the discussion more concrete, let us consider 
a well-known example of a linear quantum stochastic system in quantum optics: an optical 
cavity (see section 16.1.11 for further details of this device), shown in Figure O The cavity 
depicted in the picture is known as a standing wave or Fabry-Perot cavity and consists of one 
fully reflecting mirror at the cavity resonance frequency and one partially transmitting mirror. 
Light that is trapped inside the cavity forms a standing wave with an oscillation frequency of 
^cav, while parts of it leak through the partially transmitting mirror. The loss of light through 
this mirror is modeled as an interaction between the cavity with an incoming bosonic field A(t) 
in the vacuum state (i.e., a field with zero photons or a zero-point field) incident on the mirror. 
The dynamics for a cavity is linear and given by 

dx(t) = ~x(t)dt - y/jdA(t), 
dy(t) = yfyx(t)dtdt + dA(t), 

where 7 > is the coupling coefficient of the mirror, x(t) = (q(t),p(t)) T are the interaction 
picture position and momentum operators of the standing wave inside the cavity, and y(t) is 
the outgoing bosonic field that leaks out of the cavity. A crucial point to notice about (|2|) 
is that it is in a similar form to the classical deterministic state space representation such as 
given in ([I]), with the critical exception that ([2]) is a (quantum) stochastic system (due to the 
quantum statistical interpretation of A(t)) and involves quantities which are operator- valued 
rather than real/complex- valued. Furthermore, the system matrices A,B,C,D in ([2]) cannot 
take on arbitrary values for (|2|) to represent the dynamics of a physically meaningful system 
(see [12] and [13} Chapter 7] for further details). For instance, for arbitrary choices of A, B, C, D 
the many degrees of freedom CCR may not be satisfied for all t > as required by quantum 
mechanics; hence these matrices cannot represent a physically feasible system. In [12} I13j . 
the notion of physically realizable linear quantum stochastic systems has been introduced that 
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corresponds to open quantum harmonic oscillators (hence are physically meaningful), which 
do not include scattering processes among the bosonic fields. In particular, necessary and 
sufficient conditions have been derived on the matrices A,B,C,D for a system of the form (j2J) 
to be physically realizable. More generally, however, are linear quantum stochastic systems 
that are completely described and parameterized by three (operator-valued) parameters: its 
Hamiltonian H = \x T Rx (R £ R nxn , R = R T ), its linear coupling operator to the external 
bosonic fields L = Kx (K £ C mxn ), and its unitary scattering matrix S £ c mxm . i n particular, 
when there is no scattering involved (S = I), then it has been shown in [12] that (S,L,H) 
can be recovered from (A, B,C, D) (since S = I, here necessarily D = I) and vice-versa. 
Although [12] does not consider the scattering processes, the methods and results therein can 
be adapted accordingly to account for these processes (this is developed in section |4] of this 
paper). 

The works \12\ [T3] were motivated by the problem of the design of robust fully quantum 
controllers and left open the question of how to systematically build arbitrary linear quantum 
stochastic controllers as a suitable network of basic quantum devices. This paper addresses this 
open problem by developing synthesis results for general linear quantum stochastic systems for 
applications that are anticipated to extend beyond fully quantum controller synthesis, and it 
also proposes how to implement the synthesis in quantum optics. The organization of the rest 
of this paper is as follows. Section [2] details the mathematical modeling of linear dynamical 
quantum stochastic systems and defines the notion of an open oscillator and a generalized open 
oscillator, section [3] gives an overview of the notions of the concatenation and series product 
for generalized open oscillators as well as the concept of a reducible quantum network with 
respect to the series product, and section d] discusses the bijective correspondence between two 
descriptions of a linear dynamical quantum stochastic system. This is then followed by section [5] 
that develops the main synthesis theorem which shows how to decompose an arbitrarily complex 
linear dynamical quantum stochastic system as an interconnection of simpler one degree of 
freedom generalized open oscillators, section [6] that proposes the physical implementation of 
arbitrary one degree of freedom generalized open oscillators and direct interaction Hamiltonians 
between these oscillators, and section [7] that provides an explicit example of the application 
of the main synthesis theorem to the construction of a two degrees of freedom open oscillator. 
Finally, section [8] provides a summary of the contributions of the paper and conclusions. 

2 Mathematical modeling of linear dynamical quantum stochas- 
tic systems 

In the previous works |1CH [T2] linear dynamical quantum stochastic systems were essentially 
considered as open quantum harmonic oscillators. Here we shall consider a more general class 
of linear dynamical quantum stochastic systems consisting of the cascade of a static passive 
linear quantum network with an open quantum harmonic oscillator. However, in this paper we 
restrict our attention to synthesis of linear systems with purely quantum dynamics, whereas 
the earlier work [12] considers a more general scenario where a mixture of both quantum and 
classical dynamics is allowed (via the concept of an augmentation of a quantum linear stochastic 
system). The class of mixed classical and quantum controllers will be considered in a separate 
work. To this end, let us first recall the definition of an open quantum harmonic oscillator (for 
further details, see [T0| \12\ 115] ) ■ 

In this paper we shall use the following notations: i = \J — 1, * will denote the adjoint of a 
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linear operator as well as the conjugate of a complex number, if A = [ajk] is a matrix of linear 
operators or complex numbers, then A# = [a* k ], and A* is defined as A> = (A#) T , where T 
denotes matrix transposition. We also define = (A + A^)/2 and 9{^4} = (A — A#)/2i 

and denote the identity matrix by / whenever its size can be inferred from context and use 
Inxn to denote an n x n identity matrix. 

Let qi,pi,q2,P2, ■ ■ ■ ,qn,Pn be the canonical position and momentum operators, satisfying 
the canonical commutation relations [qj,Pk] = 2^-^, [qj,qk] = 0, \pj,Pk] = of a quantum 
harmonic oscillator with a quadratic Hamiltonian H = ^x^Rxq (xq = (q\,p\, . . . ,q n ,p n ) T ), 
where R = R T G R. The integer n will be referred to as the degrees of freedom of the oscillator. 
Digressing briefly from the main theme of this section, let us first discuss why the matrix R can 
be taken to be real (symmetric) rather than complex (Hermitian). Consider a general quadratic 
Hamiltonian H of the form H = \ Yl]=i( a jQj + PjQjPj +ljPjQj + e jP 2 j) + J2"jZl Ylk=j+i K jkqjPk, 
with ctj, €j, G R, /3j,7j G C, and [3* = V j, k, since H must be self-adjoint. Using the 
commutation relations for the canonical operators, we can then write H = ^ X^j=i( a j9j + 
^{Pj}(ljPj +PjQj) + e jPj) + YJj=i YJk=j+i KjkqjPk + c = \xqRxq + c for some real symmetric 
matrix R and a real number c = — 2 X^/=i ^{/%}- Since c contributes only a phase factor e tc 
that has no effect on the dynamics of the oscillator, as e l ^ H ~ c ^ t XQe~ l ^ H ~ c ^ t — e tHt XQe~ lHt \/t > 0, 
we may as well just discard c and take H to he H = ^x^Rxq (i.e., the original H without 
the constant term). Returning to the main discussion, let 771, . . . ,rj m be independent vacuum 
quantum white noise processes satisfying the commutation relations [rjj(t),rik(f)*] = 5jk5(t — t') 
and [rjj(t), %(i')] = Vj, k and Vi,i' > 0, and define Aj(t) = fQrjj(s)ds (j = 1, ...,m) to be 
vacuum bosonic fields satisfying the quantum Ito multiplication rules \22\ [23] 



dAj-^dA^t) = 6 jk dt, dA*(t)dA k (t) = 0, dAj(t)dA fc (t) = 0, dA*(t)dA* k {t) = 0, 

with all other remaining second order products between dAj, dA* k and dt vanishing. An open 
quantum harmonic oscillator, or simply an open oscillator, is defined as a quantum harmonic 
oscillator coupled to A(t) via the formal time- varying idealized interaction Hamiltonian [21[ 
Chapter 11] 

H Int (t) = i(L T r)(ty -tfr)(t)), (3) 

where L is a linear coupling operator given by L = Kxq with K G C mxn and r] = (rfx, ... , rj m ) T . 
Although the Hamiltonian is formal since the 77^ (i) 's are singular quantum white noise processes, 
it can be given a rigorous interpretation in terms of Markov limits (e.g., [26], [21] Chapter 11]). 
The evolution of the open oscillator is then governed by the unitary process {U (i)}t>o satisfying 
the QSDE [TO1 EH ED [27] 

dU(t) = (-iHdt + dA{t) j L - L j dA(t) - ^L j Ld?j U(t); U(0) = I. (4) 



The time-evolved canonical operators are given by x(t) = U{t)*xoU{t) and satisfy the QSDE 

dx(t) = 2e{R + S{K^K})x{t)dt + 2iO{ -ift R T ] 
x(0) = x , 

where is a canonical commutation matrix of the form = diag(J, J, . . . , J), with 



dA(t) 
dA{t)# 
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Figure 4: A generalized open oscillator. 

while the output bosonic fields Y(t) = (Yi(t), . . . ,Y n (t)) T that result from interaction of A(t) 
with the harmonic oscillator are given by Y(t) = U(t)*A(t)U(t) and satisfy the QSDE 



dY(t) 



Kx{t)dt + dA(t). 



(5) 



Note that the dynamics of x(t) and Y(t) are linear. 

The input A(t) of an open oscillator can first be passed through a static passive linear 
(quantum) network (for example, a static passive linear optical network. See, e.g., [20|, I28| for 
details) without affecting the linearity of the overall system dynamics; this is shown in Figure 
HI Such an operation effects the transformation A(t) i—> A(t) = SA(t), where S S c mxm j s a 
complex unitary matrix (i.e., S*S = SS^ = I). Thus A(t) will be a new set of vacuum noise 
fields satisfying the same Ito rule as A(t). 

Letting S = [Sjk]j,k=i,...,mi it can be shown by straightforward calculations using the quan- 
tum Ito stochastic calculus that the cascade is equivalent (in the sense that it produces the same 
dynamics for x(t) and the output y(t) of the system) to a linear quantum system whose dynam- 
ics is governed by a unitary process {U(t)}t>o satisfying the QSDE (for a general treatment, 
see [29]) 



dU{t) = ( J2{S jk -5 jk )dA jk {t)-iHdt + dA(tfL-L j SdA(t) 



(6) 



-L ] Ldt ] U{t);U(0) =/, 



where Aj k (t) (j, k = 1, . . . , m) are fundamental processes, called the gauge processes, satisfying 
the quantum Ito rules 

dk jk (t)dkj> k >(t) = 5 kf dA jk ,(t), dAj(t)dA kl (t) = S jk dAi(t), dA jk dAi(t)* = 5 kl dA*{t), 

with all other remaining second order products between dAji(t) and dAji (t) , dAf, (t) , dt van- 
ishing. This yields the following dynamics for x(t) = U(t)*xoU(t) and the system output 
y(t) = U(t)*A(t)U(t): 



with 



dx(t) 
dy(t) 

A 
B 
C 
D 



Ax(t)dt + B 



dA(t) 
dA{t)# 

Cx{t)dt + DdA(t), 



2Q(R + S{K j K}), 
2iS[ -K^S K T S# 
SK, 
S. 



(7) 
(8) 

(9) 
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For convenience, in the remainder of the paper we shall refer to the cascade of a static passive 
linear quantum network with an open oscillator as a generalized open oscillator. 

Let G be a generalized open oscillator that evolves according to the QSDE ([6]) with given 
parameters S, L = Kxo, and H = ^x^Rxq. For compactness, we shall use a shorthand 
notation of and denote such a generalized open oscillator by G = (S,L,H). In the next 
section we briefly recall the concatenation and series product developed in [29] that allows 
one to systematically obtain the parameters of a generalized open oscillator built up from an 
interconnection of generalized open oscillators of one degree of freedom. 

3 The concatenation and series product of generalized open 
oscillators and reducible quantum networks 

In this section we will recall the formalisms of concatenation product, series product, and 
reducible networks (with respect to the series product) developed in [29\ for the manipulation 
of networks of generalized open oscillators as well as more general Markov open quantum 
systems. 

Let G\ = (Si, Kixifi, ^xf RiXifi) and G2 = (5*2, K2X 2 ,o, \x^ R 2 X2fi) be two generalized 
open oscillators, where Xf~Q = Xfc(0). The concatenation product G\ EE G2 of G\ and G2 is 
defined as 

(rrt X r f 1 on \ 

Sim, [Kixifi, K 2 X2fl) , -x lfi Rix lfi + -x 2)0 -R2^2,o J , 



where 



S\m2 



St 
S 2 



It is important to note here that the possibility that x± = ^2,0 or that some components of 
xi t o coincide with those of 212,0 are allowed. If Gi and G2 are independent oscillators (i.e., the 
components of x^o act on a distinct Hilbert space to that of the components of x 2j o), then the 
concatenation can be interpreted simply as the "stacking" or grouping of the variables of two 
noninteracting generalized open oscillators to form a larger generalized open oscillator. 

It is also possible to feed the output of a system G\ to the input of system G2, with the 
proviso that Gi and G 2 have the same number of input and output channels. This operation 
of cascading or loading of G2 onto G\ is represented by the series product G2 < G\ defined by 

G 2 <Gi = fs 2 Si, K 2 x 2 fl + S^A'ix^o, ^x\ fl Rix lfi 

+ ^2,0^2,0 + Yi X l°( K * S2Kl ~ KlS*K*) Xl ^ . 

Note that G2 < G\ is again a generalized open oscillator with a scattering matrix, coupling 
operator, and Hamiltonian as given by the above formula. 

With concatenation and series products having been defined, we now come to the important 
notion of a reducible network with respect to the series product (which we shall henceforth refer 
to more simply as just a reducible network) of generalized open oscillators. This network 
consists of I generalized open oscillators Gk = (Sk, Lk, Hk), with Lf- = K^xtfl and H^ = 
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ixj^Q-RfcXj^o, k = along with the specification of a direct interaction Hamiltonian 

H d = J2j ^2k=j+i x JfiRjk x k,o {Rjk G M 2x2 ) and a list S = {Gk<Gj} of series connections among 
generalized open oscillators Gj and Gk, j ^ k, with the condition that each input and each 
output has at most one connection, i.e., lists of connections such as {G2 <G%, G3 <G2, G\ <G^} 
are disallowed. Such a reducible network Af again forms a generalized open oscillator and 
is denoted by Af = {{Gk}k=i,...,i, H d , S}. Note that if Afo is a reducible network defined 
as Ao = {{Gk}k=i,...,h 0, S} = (Sq, Lq, Hq), then Af, which is Afo equipped with the direct 
interaction Hamiltonian H d , is simply given by Af = Afo EH (0, 0, H) = (So, Lq, Hq + H d ). 

The notion of a reducible network was introduced in [29] to study networks that are free 
of "algebraic loops" such as when connections like {G2 < Gi,G% <\G2,G\ < G3} are present. 
The theory in [29] is not sufficiently general to treat networks with algebraic loops; they can 
be treated in the more general framework of quantum feedback networks developed in [30j . 
Since this work is based on |29j, we also restrict our attention to reducible networks, but as 
we shall show in section [5] this is actually sufficient to develop a network synthesis theory of 
linear quantum stochastic systems. A network synthesis theory can indeed also be developed 
using the theory of quantum feedback networks of [30] , and this has been pursued in a separate 
work [3Tj . 

Two important decompositions of a generalized open oscillator based on the series product 
that will be exploited in this paper are 

(S,L,H) = (I, L, H) < (S, 0, 0), (10) 
(S,L,H) = (S,0,0)«(I,^L,H), (11) 

where (S, 0, 0) represents a static passive linear network implementing the unitary matrix S. 

4 Correspondence between system matrices (A, B, C, D) and 
the parameters S, L, H 

In |12] it has been shown that for S = I, then D = I, and there is a bijective correspondence 
between the system matrices (A, B,C, I) of a physically realizable linear quantum stochastic 
system [12\ section III] and the parameters K, R of an open oscillator; see Theorem 3.4 therein 
(however, note that the B, C, and D matrices are defined slightly differently from here because 
[12j expresses all equations in terms of quadratures of the bosonic fields rather than their 
modes). Here we shall show that allowing for an arbitrary complex unitary scattering matrix 
S, a bijective correspondence between the system parameters (A, B,C, D) of an extended 
notion of a physically realizable linear quantum stochastic system and the parameters S, K, R 
of a generalized open oscillator (in particular, D = S) can be established. We begin by noting 
that we may write the dynamics ([8]) in the following way: 

y(t) = Sy'(t), 

with y'(t) defined as 

d y '(t) = S ] Kx(t)dt + dA(t). 

Then by defining K 1 = S^K and substituting K = SK' in ([9|), we see that x(t) in ([7]), and 
y'(t) constitutes the dynamics for the open oscillator (I,K ! xq, ^x^Rxo) with system matrices 
given by (A,B,S^C,I). Since D = S and (S,L,H) = (S,0,0) < (I,K'x , ^Rx ) (cf. ([II])), 
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from [121 Theorem 3.4] we see that there is a bijective correspondence between (A, B, S^C) and 
the parameters (K',R) and that one set of parameters may be uniquely recovered from the 
other. Therefore, we may define a system of the form ([2]) to be physically realizable (extending 
the notion in [12] ) if it represents the dynamics of a generalized open oscillator (this idea 
already appears in [131 Chapter 7]; see Remark 7.3.8 therein). This implies that a system ([2]) 
with matrices (^4, B, C, D) is physically realizable if and only if D is a complex unitary matrix 
and (A, B, D^C, I) are the system matrices of a physically realizable system in the sense of |12| 
(i.e., (A, B, D^C, I) are the system matrices of an open oscillator). Therefore, we may state 
the following theorem. 

Theorem 1 There is a bijective correspondence between the system matrices (A, B, C, D) and 
the parameters (S,K,R) of a generalized open oscillator. For given (S,K,R), the correspond- 
ing system matrices are uniquely given by ([9]). Conversely, for given (A,B,C,D), which are 
the system matrices of a generalized open oscillator G with parameters (S,K,R), then D 
is unitary, and S = D and (A,B,D^C,I) are the system matrices of some open oscillator 
G' = (I,K'xq,^XqRxq). The parameters (K',R) of the open oscillator G' is uniquely deter- 
mined from (A,B,D^C) by \12\ Theorem 3.4] (by suitably adapting the matrices B and D^C), 
from which the parameter K of G is then uniquely determined as K = DK' . 

Due to this interchangeability of the description by (^4, B, C, D) and by (S, K, R) for a 
generalized open oscillator, it does not matter with which set of parameters one works with. 
However, for convenience of analysis in the remainder of the paper we shall work exclusively 
with the parameters (S,K,R). 

5 Main synthesis theorem 

Suppose that there are two independent generalized open oscillators coupled to m independent 
bosonic fields, with m output channels: an n\ degrees of freedom oscillator G\ = (Si, L\, H\) 
with canonical operators x\ = (qi,i,Pi,i, • • • , qi,n 1 ,Pi,n 1 ) T , Hamiltonian operator Hi = ^x T R\x±, 
coupling operator L\ = K\x\, and scattering matrix Si, and, similarly, an n2 degrees of free- 
dom oscillator G2 = (6*2,^2,-^2) with canonical operators X2 = (q2,i,P2,i, ■ ■ ■ , l2,n 2 iP2,n 2 ) T ^ 
Hamiltonian operator H2 = \x\R2X2, coupling operator L2 = -^2^2) and unitary scattering 
matrix S2. 

Consider now a reducible quantum network A/12 constructed from G\ and G2 as A/12 = 
{{Gi, G2}, Hf 2 , G2 <G±}, as shown in Figure [5j where Hf 2 is a direct interaction Hamiltonian 
term between G\ and G2 given by 

Hf 2 = ^xjR 12 x 2 + \xlR^ 2Xl - \jAS2U - L\S\L 2 ); R n G M 2x2 
= x T 2 R\ 2 x x - \jAS2Lx - L T 2 S*L*) 
= x^Rj 2 -^- i (KlS2K 1 -K 2 T S*K*) ) jx 1 , 

where we recall that A& denotes the elementwise adjoint of a matrix of operators A and the 
second equality holds, since elements of L\ commute with those L 2 . Also note that the matrix 
2 ^(K 2 S2Ki — K 2 ' S 2 Kf") is real. Some straightforward calculations (see [29] for details) then 
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show that we may write 
M 2 



{S 2 S U S 2 L Y + L 2 , Hi + H 2 + H[ 2 + Hf 2 ), 



where H( 2 = 2 ^(L 2 S 2 Li — L\S 2 L 2 ). Now let us look closely at the Hamiltonian term of A/"i 2 . 
Note that after plugging in the definition of H\, H 2 , Hf 2 , and H( 2 , we may write 



H\ + H 2 + h( 2 + Hf 2 



T 



i?l 
#12 



-^12 
#2 



,ri 
x 2 



Letting x = {x[ ,x 2 ) , SWi 



S^Si, and defining 
Ri 



R 
K 



R 12 



we see that 



A/i 



12 



-R12 
#2 

: [ S 2 #l *T 2 ], 
1 T 



(12) 
(13) 

(14) 



Kx, and-ff^-i 



Therefore, Af 12 = (S 2< -i, L 2< -i, H 2 <-i), with 5 2< _i = S^Si, L 2 ^ 
In other words, a reducible network formed by a bilinear direct interaction and cascade con- 
nection of two generalized open oscillators having the same number of input and output fields 
results in another generalized open oscillator with a degrees of freedom which is the sum of the 
degrees of freedom of the two constituent oscillators and having the same number of inputs 
and outputs. 

By repeated application of the above construction, we can prove the following synthesis 
theorem. 



Theorem 2 Let G be an n degrees of freedom generalized open oscillator with Hamiltonian 
matrix R £ ]j2nx2n^ coupling matrix K £ ^rnx2n^ an g un n ar y scattering matrix S £ C mxm . 
Let R be written in terms of blocks of 2x2 matrices as R = [Rjklj k=i ... m where the Rjk 's are 
real 2x2 matrices satisfying Rkj = Rj k for all j, k, and let K be written as 

K=[K l K 2 ... K n ], 

where, for each j, Kj £ C mx2 . For j = 1, . . . , n, let Gj = (Sj, KjXj, ^xjRjjXj) be independent 
one degree of freedom generalized open oscillators with canonical operators Xj = (qj,Pj) T , m 
output fields, Hamiltonian matrix Rjj, coupling matrix Kj, and scattering matrix Sj. Also, 




Figure 5: Cascade connection of G± and G 2 with indirect interaction Hf, 
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define S k ^j for j < k + 1 as S k ^j = T\* =j Si = S k --- S j+1 Sj for j < k, S k< ^ k = S k , and 
Sfc^+i = Imxm, and let H d be a direct interaction Hamiltonian given by 

n— 1 n / 1 \ 

Hd = E E s * fl l - - ^t i+ i^f) (15) 

j=l k=j+l v J 

If Si,...,S n satisfies S n S n _i ■ ■ ■ Si = S and K k satisfies K k = Sl^ k+1 K k for k = 1, . . . , n, 
then the reducible network of harmonic oscillators N given by N = {{Gi, . . . , G n }, H d , {G2 < 
G±, G3 < G2, ■ ■ ■ , G n <\ G n _i}} is equivalent to G. That is, G can be synthesized via a series 
connection G n <\. . .<G2<G\ of n one degree of freedom generalized open oscillators, along with 
a suitable bilinear direct interaction Hamiltonian involving the canonical operators of these 
oscillators. In particular, if S = I m xm (no scattering), then S k can be chosen to be S k = I m xm 
and K k can be chosen to be K k = K k for k = 1, . . . , n. 

Proof. Let Hj = RjjXj, Lj = KjXj and 

H l = E f L ) E - E A^rih] > k * 2 - 

3=2 \ 1=1 1=1 J 

Let us begin with the series connection G12 = G2 < G\. By analogous calculations as given 
above for the two oscillator case, it is given by 

G12 = (S2S1, S2L1 + L2, H\ + H2 + i?2 ). 

Repeating this calculation recursively for G123 = G3 < G12, G1234 = G4 < G123,..., 
Gi2... n = G n < G 12 ...( n _i), we obtain at the end that 

/ n n \ 

Gi2... n = I £n«-i>E S n «_fe + iLjfc, H k + H[ J . 



k=l k=l 



Noting that Hn may be rewritten as 



1 n—l n 

Hf n = ^ E E (4wi^ - 
j=i fc=j+i 

.. n— 1 n 

= ^ E E - L l s U + i L f) 

j=i k=j+i 

1 n— 1 n 

= ^E E ^KlSk^+iKj-K^St^KpXj, 
j=i k=j+i 

where the second equality holds since Lj commutes with L k whenever j 7^ k, we find that 

n j n n ^ 

J2 H k + Hl + H d = -J2J2xjR jk x k = -x T Rx, x = (x[,xl,...,x T n ) T ■ 

k=l j=l k=l 
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Therefore, if S\, . . . , S n and K\, . . . , K n satisfy the conditions stated in the theorem, we find 
that TV = {{G 1} . . . , G n }, H d , {G 2 < G 1} G 3 < G 2 , . . . , G n < G„_i} is given by 



That is, /V is a linear quantum stochastic system with Hamiltonian matrix R, coupling matrix 
K, and scattering matrix S, and is therefore equivalent to G. This completes the proof of the 



Therefore, according to the theorem, synthesis of an arbitrary n degrees of freedom linear 
dynamical quantum stochastic system is in principle possible if the following two requirements 
can be met: 

1. Arbitrary one degree of freedom open oscillators G = (I, L, H) with m input and output 
fields can be synthesized. In particular, it follows from this that one degree of freedom 
generalized open oscillators G = (S, L, H) can be synthesized as G = (I, L, H ) < (S, 0, 0). 

2. The bilinear interaction Hamiltonian H d as given by (|15p can be synthesized. 

One can observe certain parallels between the quantum synthesis described in the theorem 
with the active state-space synthesis method of linear electrical network synthesis theory (e.g., 
[191 Chapter 13]). To begin with, we may think of each oscillator Gj as a type of noisy 
quantum integrator, as the counterpart of a classical integrator (built from an operational 
amplifier, resistors, and capacitors) in an electrical network, and in both settings synthesis can 
be achieved by suitably cascading these basic integrating components. We may also view the 
direct interaction Hamiltonian between two oscillators as acting like a type of mutual "state 
feedback" between the oscillators, much like the state feedback employed in electrical network 
synthesis. However, because of the quite distinct nature of classical and quantum systems, 
of course the parallels should not be taken to be "exact" or "precise" in any way, the nature 
of these parallels are in spirit rather than detail. Certainly, classical active synthesis theory 
cannot be applied directly to linear quantum stochastic systems because of certain constraints 
that a noisy quantum integrator must satisfy that are not required of its classical counterpart, 
and the classical theory is deterministic while in the quantum theory, quantum stochastic 
noise plays a crucial role, for instance, to preserve the canonical commutation relations in 
open quantum systems. To highlight another significant difference between the two physical 
systems, we note that losses in linear electrical systems may be modeled by inserting resistors 
as dissipative components of the system, while in linear quantum systems, losses are modeled 
by lossy couplings to quantum noises (heat baths). 

6 Systematic synthesis of linear quantum stochastic systems 

This section details the construction of arbitrary one degree of freedom open quantum harmonic 
oscillators and implementation of bilinear direct interactions among the canonical operators of 
these oscillators, at least approximately, in the context of quantum optics, using various linear 
and nonlinear quantum optical components. 

We begin with a description of some key quantum optical components that will be required 
for the synthesis. This is followed by a discussion of general synthesis of one degree of freedom 
open oscillators and finally by a discussion of the implementation of bilinear direct interaction 
Hamiltonians among different one degree of freedom open oscillators. 




synthesis theorem. 
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Figure 6: Two types of optical cavities: a standing wave or Fabry-Perot cavity (left) and a 
(three mirror) ring cavity (right). Arrows indicate the direction of propagation of light in the 
cavity. Black rectangles denote fully reflecting mirrors, while white rectangles denote partially 
transmitting mirrors. 



6.1 Essential quantum optical components 

6.1.1 Optical cavities 

An optical cavity is a system of fully reflecting or partially transmitting mirrors in which a 
light beam is trapped and is either bounced repeatedly from the mirrors to form a standing 
wave or circulates inside the cavity (as in a ring cavity); see Figure If there are transmitting 
mirrors present, then light can escape or leak out from the cavity, introducing losses to the 
cavity. 

A cavity is mathematically modeled by a Hamiltonian H cav = oj cav a*a, where oj cav is the 
resonance frequency of the cavity and a = is the (non-self-adjoint) cavity annihilation 
operator or cavity mode satisfying the commutation relation [a, a*] = 1. Here q = a + a* is the 
position operator of the cavity mode (also called the amplitude quadrature of the mode) and 
p = —ia + ia* is the momentum operator of the cavity mode (also called the phase quadrature 
of the mode). If there is a transmission mirror, say, M, then losses through this mirror are 
modeled as having a vacuum bosonic noise field A(t) incident at this mirror and interacting 
with the cavity mode via the idealized Hamiltonian Hj n t given in ([3]) with L = \fna, where k is 
a positive constant called the mirror coupling coefficient. When there are several leaky mirrors, 
then the losses are modeled by a sum of such interaction Hamiltonians, one for each mirror and 
with each mirror having its own distinct vacuum bosonic field. The total Hamiltonian of the 
cavity is then just the sum of H cav and the interaction Hamiltonians. More generally, the field 
incident at a transmitting mirror need not be a vacuum field, but can be other types of fields, 
such as a coherent laser field. Nonetheless, the interaction of the cavity mode with such fields 
via the mirror will still be governed by (j3J) with a coupling operator of the form L = ^JUa. 

6.1.2 Degenerate parametric amplifier 

In order to amplify a quadrature of the cavity mode, for example, to counter losses in that 
quadrature caused by light escaping through a transmitting mirror, one can employ a non- 
linear optical crystal and a classical pump beam in the configuration of a degenerate parametric 
amplifier (DPA), following the treatment in \21\ section 10.2]. The pump beam acts as a source 
of additional quanta for amplification and, in the nonlinear crystal, an interaction takes place 
in which photons of the pump beam are annihilated to create photons of the cavity mode. In 
an optical cavity, such as a ring cavity shown in Figure we place the crystal in one arm of 
the cavity (for example, in the arm between mirrors Ml and M2) and shine the crystal with 
a strong coherent pump beam of (angular) frequency u p given by uj p = 2oj r , where uj r is some 
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Classical pump Classical pump 




Ml 



Figure 7: A DPA consisting of a classically pumped nonlinear crystal in a three mirror ring 
cavity. 



Figure 8: Schematic representation of a DPA. The white rectangle symbolizes the nonlinear 
crystal, while the diagonal arrow into the rectangle denotes the pump beam. 

reference frequency. Here the mirrors at the end the arms should be chosen such that they 
do not reflect light beams of frequency u p . A schematic representation of a DPA (a nonlinear 
crystal with a classical pump) is shown in Figure [8l 

Remark 3 In the remaining figures, black rectangles will be used to denote mirrors which are 
fully reflecting at the cavity frequency and fully transmitting at the pump frequency (whenever 
a pump beam is employed) , while white rectangles denote partially transmitting mirrors at the 
cavity frequency. 

Let a = ^4^2 be the cavity mode, and let the cavity frequency <jj cav be detuned from oj r and 
given by u} cav = u; r + A, where A is the frequency detuning. The crystal facilitates an energy 
exchange interaction between the cavity mode and pump beam. By the assumption that the 
pump beam is intense and not depleted in this interaction, it may be assumed to be classical, 
in which case the crystal-pump-cavity interaction can be modeled using the (time-varying) 
Hamiltonian H(t) = oj cav a*a + |(ee-^*(a*) 2 - e*e i ^ t a 2 ) [2U equation 10.2.1], where e is a 
complex number representing the effective pump intensity. By transforming to a rotating frame 
with respect to w r = y (i-e., by application of the transformation a i— > ae l ~^ 1 ; see |21| section 
10.2.1] for a derivation of the equations of motion of the DPA in the rotating frame), H can 
be reexpressed as 




e a ) 



and be written compactly as H = ^XqRxq + c (recall xq = (q,p) T ), where 



2 I ^ + A-i(e-0 



(16) 



and c is a real number. Since c merely contributes a phase factor that has no effect on the 
overall dynamics of the system operators, it plays no essential role and can simply be ignored 
(cf. section [2]). Note that transformation to a rotating frame effects the following: If a(t) is the 
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Figure 9: Two cavity modes interacting in a single classically pumped nonlinear crystal. The 
diagonal arrow into the crystal denotes the pump beam. 




Figure 10: Schematic representation of a two-mode squeezer. 



evolution of a under the original time-varying Hamiltonian H(t) = Lu cav a*a + |(ee l ^p*( a *) 2 — 
e*e lulpt a 2 ) and we define a(t) = a(t)e* avt (i.e., a(t) is a(t) in a frame rotating at frequency 
u r ), then a(t) coincides with the time evolution of a under the time-independent Hamiltonian 
H = t^XqRxq. In other words, in this rotating frame, the DPA can be viewed as a harmonic 
oscillator with quadratic Hamiltonian ^XqRxq- 



6.1.3 Two-mode squeezing 

If two cavities are positioned in such a way that the beams circulating in them intersect one 
another, then these beams will merely pass through each other without interacting. One way 
of making the beams interact is to have their paths intersect inside a x nonlinear optical 
crystal. Typically, to facilitate such an interaction, one or two auxiliary pump beams are also 
employed as a source of quanta/energy. For instance, in a x^ optical crystal in which the 
modes of two cavities interact with an undepleted classical pump beam as depicted in Figure 
[9l the interaction can be modeled by the Hamiltonian 

H{t) = ^(ee'^alal - eV^a^), 

where e is a complex number representing the effective intensity of the pump beam and oj p is the 
pump frequency. Transforming to a rotating frame at half the pump frequency by applying the 

■ ^p , ■ 

rotating frame transformation a± i— > a\e l ^~ and a,2 i— >■ a^e*" , H(t) can be expressed in this 
new frame in the time-invariant form H = ^(ea^a^ - e * a i a 2)- This type of Hamiltonian is called 
a two-mode squeezing Hamiltonian, as it simultaneously affects squeezing in one quadrature of 
(possibly rotated versions of) a\ and and will play an important role later on in the paper. 
A two-mode squeezer is schematically represented by the symbol shown in Figure [TU1 

Remark 4 It will be implicitly assumed in this paper that the equations for the dynamics of 
generalized open operators are given with respect to a common rotating frame of frequency 
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u r , including the transformation of all bosonic noises Ai(t) according to Ai(t) h- > Ai(t)e luJrt , 
and that classical pumps employed are all of frequency uj v = 2u r . This is a natural setting in 
quantum optics where a rotating frame is essential for obtaining linear time invariant QSDE 
models for active devices that require an external source of quanta. In a control setting, this 
means both the quantum plant and the controller equations have been expressed in the same 
rotating frame. 



6.2 Static linear optical devices and networks 

Static linear optical devices implement static linear transformations (meaning that the trans- 
formation can be represented by a complex square matrix) of a set of independent incoming 
single mode fields, such as the field in a cavity, a = (a%, 02, . . . a m ) T to an equal number 
a' = (dj, a' 2 , . . . a' m ) T of independent outgoing fields. The incoming fields satisfy the commuta- 
tion relations [aj,a k ] = and [aj,a£] = 5j k . The incoming fields may also be vacuum bosonic 
fields A(t) = (A\(t),A2(t), . . . , A m {t)) T with outgoing bosonic fields (that need no longer be in 
the vacuum state) A'(t) = (A\(t), A' 2 (t), . . . , A' m (t)) T . In the latter, the commutation relations 
are [dAj(t),dA k (t)] = and [dAj(t), dAf.(t)*] = Sji-dt. However, to avoid cumbersome and 
unnecessary repetitions, in the following we shall only discuss the operation of a static linear 
optical device in the context of single mode fields. The operation is completely analogous for 
bosonic incoming and outgoing fields and requires only making substitutions such as a — > A(t), 
a' - 
etc. 



A'(t), [aj,a k ] = -)• [dAj(t),dA k (t)] = 0, and 



5 jk -> [dAj(t),dA k (t)*] = 8 jk dt, 



The operation of a static linear optical device can mathematically be expressed as 



7 '# 



Q 



-,# 



; Q 



Q2 

Q2 Q\ 



Qi 

r 2 



where Qi,Q2 £ C mxm and S is a quasi-unitary matrix |201 section 3.1] satisfying 



Q 



1 
-/ 



1 
-/ 



A consequence of the quasi-unitarity of Q is that it preserves the commutation relations among 
the fields, that is, to say that the output fields a' satisfy the same commutation relations as 
a. Another important property of a quasi-unitary matrix is that it has an inverse Q~ l given 
by Q' 1 = GQ^G, where G = J _° 7 , and this inverse is again quasi-unitary, i.e., the set of 
quasi- unitary matrices of the same dimension form a group. 

In the case where the submatrix Q2 of Q is Q2 = 0, the device does not mix creation and 
annihilation operators of the fields, and it necessarily follows that Q\ is a complex unitary 
matrix. Such devices are said to be static passive linear optical devices because they do not 
require any external source of quanta for their operation. It is well known that any passive 
network can be constructed using only beam splitters and mirrors (e.g., see references 2-4 
in |28j). In all other cases, the devices are static active. Specific passive and static devices that 
will be utilized in this paper will be discussed in the following. 



6.2.1 Phase shifter 

A phase shifter is a device that produces an outgoing field that is a phase shifted version of 
the incoming field. That is, if there is one input field a, then the output field is a' = e %e a 
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Figure 11: Phase shifter with a phase shift of 9 radians. 



for some real number 9, called the phase shift; a phase shifter is schematically represented by 
the symbol shown in Figure [TTJ By definition, a phase shifter is a static passive device. The 
transformation matrix Qps of a phase shifter with a single input field is given by 



Q 



o 



PS 



6.2.2 Beam splitter 



A beam splitter is a static and passive device that forms a linear combination of two input fields 
a\ and 02 to produce two output fields a'i and a' 2 such that energy is conserved: a\ai + a\a2 = 
(a^Ya'i + (a 2 )*a 2 . The transformation affected by a beam splitter can be written as 



Qbs 



B 
B# 



where B is a unitary matrix given by 
B = e £ / 2 



e-^/ 2 



cos(G)/2 sin(G)/2 
-sin(G)/2 cos(9)/2 











Here H,0,<1?,\I/ are real numbers. is called the mixing angle of the beam splitter, and it is 
the most important parameter. <3? and ^ introduce a phase difference in the two incoming and 
outgoing modes, respectively, while H introduces an overall phase shift in both modes. 

A particularly useful result on the operation of a beam splitter with H = ^ = $ = 0is 
that it can be modeled by an effective Hamiltonian H BS given by H BS = i®{a\a2 — aia^) 
(see [20], section 4.1] for details). This means that in this case we have the representation 



Qbs 



-,# 



exp(iH° BS ) 



,# 



exp(-iH BS ), 



where a = (01,02) ■ More generally, it follows from this, by considering phase shifted inputs 



• s+# 



a\e l 2 and 02 



a2e* 2 (J) being an arbitrary real number), that a beam splitter with 



and ^ = — $ will have the effective Hamiltonian Hbs = i®(e l<s a\a2 — e^aiOrj) 



aa*a2 + a*aia2, with a = iQe~ l ® . This is the most general type of beam splitter that will be 
employed in the realization theory of this paper. A beam splitter with a Hamiltonian of the 
form Hbs is represented schematically using the symbol in Figure [T2l 



6.2.3 Squeezer 

Let there be a single input mode a. Write a as a = q' + ip' , where q' = (a + a*)/ 2 is the real 
or amplitude quadrature of a and p' = (a — a*)/2i is the imaginary or phase quadrature of a. 
Squeezing of a field is an operation in which the variance of one quadrature, either q' or p' , is 
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Figure 12: Schematic representation of a beam splitter. 




Figure 13: Schematic representation of a squeezer. 



squeezed or attenuated (it becomes less noisy) at the expense of increasing the variance of the 
other quadrature (it becomes noisier). A device that performs squeezing of a field is called a 
squeezer. An ideal squeezer affects the transformation Qsqueezer given by 



Q 



squeezer 



cosh(s) 



sinh(s) 



e sinh(s) cosh(s) 



where s and 9 are real parameters. We shall refer to s as the squeezing parameter and 9 as 
the phase angle. For s < 0, the squeezer squeezes the amplitude quadrature of e~ l ^a (a phase 
shifted version of a) while if s > 0, it squeezes the phase quadrature and then shifts the phase 
of the squeezed field by |. A squeezer with parameters s, 9 is schematically represented by the 
symbol shown in Figure [TBI 

A squeezer can be implemented, for instance, by using a combination of a parametric 
amplifier and a beam splitter for single mode fields \20\ section 6.1] or as a DPA with a 
transmitting mirror for bosonic fields [21} section 7.2.9]. It is easy to see that Q~^q Ueezer is given 
by 



Q 



squeezer 



cosh(— s) 
~ %e sinh(— s) 



sinh(- 



cosh(— s) 

6.2.4 Static optical linear networks 



cosh(s) 
-e~ t0 sinh(s) 



jo 



sinh(s) 



cosh(s) 



It is known that an arbitrary static linear optical network can be decomposed as a cascade of 
simpler networks. In particular, any quasi- unitary matrix Q can be constructively decomposed 
as [28]: 



Q 





" Ai 







D ' 




" ^3 





exp 


A* _ 


exp 


D 





exp 








exp A\ 






exp^4f 



cosh D 
sinh D 



sinh D 
coshD 



exp A3 




exp 



4 



where A\ and A% are skew symmetric complex matrices and D is a real diagonal matrix. The 
first and third matrix exponential represent passive static networks that can be implemented by 
beam splitters and mirrors, while the second exponential represents an independent collection 
of squeezers (with trivial phase angles) each acting on a distinct field. 

In summary, in any static linear optical network the incident fields can be thought of as 
going through a sequence of three operations: they are initially mixed by a passive network, 
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Figure 14: Implementation of a squeezer with arbitrary phase angle employing a squeezer with 
a zero phase angle and two phase shifters. 



then they undergo squeezing, and finally they are subjected to another passive transformation. 
In the special case where the entire network is passive, the squeezing parameters (i.e., elements 
of the D matrix) are zero. 

For example, a squeezer with arbitrary phase angle 6 can be constructed by sandwiching a 
squeezer with phase angle between a —9/2 phase shifter at its input and a 9/2 phase shifter 
at its output, respectively. This is shown in Figure [Til 

6.3 Synthesis of one degree of freedom open oscillators 

One degree of freedom open oscillators are completely described by a real symmetric Hamil- 
tonian matrix R = R T G K 2x2 and complex coupling matrix K G C mx2 . Thus one needs 
to be able to implement both R and K. Here we shall propose the realization of one degree 
of freedom open quantum harmonic oscillators based around a ring cavity structure, such as 
shown in Figure [6l using fully reflecting and partially reflecting mirrors and nonlinear optical 
elements appropriately placed between the mirrors. 

The matrix R determines the quadratic Hamiltonian H = ^x T Rx and in a one-dimensional 
setup such a quadratic Hamiltonian can be realized with a DPA as discussed in section 16.1.21 
From (|16p . it is easily inspected that any real symmetric matrix R can be realized by suitably 
choosing the complex effective pump intensity parameter e and the cavity detuning parameter 
A of the DPA. In fact, for any particular R, the choice of parameters is unique. For example, 
to realize 



one solves the set of equations 

A - 9{e} = 2, 3fi{e} = -4, and A + 9{e} = 1 

for A, ^{e}, 3?{e} to yield the unique solution A = 3/2 and e = — 4 — i/2. 

Now, we turn to consider realization of the coupling operator L = Kxq. Let us write 
K = [ Kj ... ] T , where K\ G C lx2 for each / = l,...,m. Each K[ represents the 

coupling of the oscillator to the bosonic field Ai, and so it suffices to study how to implement 
a single linear coupling to just one field. To this end, suppose now that there is only one field 
A(t) coupled to the oscillator via a linear coupling operator L = Kxq for some K G C lx2 . It 
will be more convenient to express L = aq + ftp in terms of the oscillator annihilation operator 
a and creation operator a* defined by a = (q + ip)/2 and a* = (q — ip)/2. Therefore, we write 
L = aa + (3a*, with a = and /3 = . Consider the physical scheme shown in Figure [151 
partly inspired by a scheme proposed by Wiseman and Milburn for quantum nondemolition 
measurement of the position operator, treated at the level of master equations |32j (whereas 
here we consider unitary models and QSDEs). In this scheme, additional mirrors are used to 
implement an auxiliary cavity mode b of the same frequency as the reference frequency co r (cf. 
Remark H]) . The auxiliary cavity b interacts with a via a cascade of a two-mode squeezer and 



LINEAR QUANTUM STOCHASTIC NETWORK SYNTHESIS 



23 



Aft) 



Yft) 



Aft) 

JL 



* 1 Ai iv i haru 
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Figure 15: Scheme for (approximate) implementation of a coupling L = da + /3a* to cavity 
mode a using an auxiliary cavity b (whose dynamics is adiabatically eliminated), a two- mode 
squeezer, and a beam splitter with the appropriate parameters. The left figure is a block 
diagram showing the fast mode b interacting with the slow mode a via the direct interaction 
Hamiltonian H a ^, implemented by the two-mode squeezer and the beam splitter, and also 
interacting with a 180° phase shifted input field A(t) to produce the output field Y(t). The 
right figure details the physical implementation of the block diagram. 



a beam splitter. The combination of the nonlinear crystal and beam splitter implements an 
overall interaction Hamiltonian H ao , in a rotating frame at frequency oj r (equal to half the 
pump frequency of the two-mode squeezer) , of the form 

H ab = ^(eia*b* - e\ab) + l -(e 2 a*b - e* 2 ab*), (17) 

where e\ is the effective pump intensity of the two- mode squeezer and e 2 is given by e 2 = 20e -i *, 
where is the mixing angle of the beam splitter and $ is the relative phase introduced between 
the input fields by the beam splitter. Assuming that the coupling coefficient 72 of the partially 
transmitting mirror M on b is such that b is heavily damped compared to a, b will have 
much faster dynamics than a and thus allows one to adiabatically eliminate b to obtain a 
reduced dynamics for a only. A rigorous foundation for such adiabatic elimination or singular 
perturbation procedure has recently been developed in [33]. Based on this theory, the adiabatic 
elimination results developed in Appendix [8] show that, after the elimination of b, the resulting 
coupling operator to a will be given by 

L = -^—(-e* 2 a + eia*). 

Therefore, it becomes clear that by choosing the parameters €1,62,72 with 72 large and such 
that 

a = — % and /3 = — =, (18) 

it is possible to approximately implement any coefficients d and (3 in a linear coupling operator 
L = da + (3a* . Note that a ir radian phase shifter in front of A in Figure [15] is required 
to compensate for the scattering term in the unitary model that is obtained after adiabatic 
elimination (cf. Appendix [8]). 

Moreover, for the special case where d,f3 satisfy d is real and d > > we also propose 
an alternative implementation of the linear coupling based on preprocessing and postprocessing 
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with squeezed bosonic fields (see Appendix [8] for details). To this end, let 7 
and consider the interaction Hamiltonian 

H Int (t) = i(L V (ty - L*n(t)) 

= i((aa + Pa*)r}{t)* - (da* + f3*a)r](t)). 



a 2 



> 0, 



Let us rewrite this Hamiltonian as follows 

Hlntit) 



i(a(arj(t)* - 



P* V (t))-a*(a V (t) 
-a*r,'(t)), 



HtT)) 



where r]'(t) = 
PA(t)*), and 



(ar](t) -f3ri(t)*). Letting Z(t) = f*r)'(s)ds, we have that Z{t) = -±={&A(t) 



Z(t) 
Z{t)* 



Q 



A(t) 
A(t)* 



Q 



a 

V7 







a 


V7 





The main idea is that instead of considering an oscillator interacting with A(t), we consider the 
same oscillator interacting with the new field Z(t) via the interaction Hamiltonian Hj n t(t) = 
iy/j(arf(t)* — a*rf(t)). Since a 2 - \/3\ 2 = 7 > 0, we see that {a/^) 2 - l/V^/rl 2 = 1> ^ vom which 
it follows that Q is a quasi-unitary linear transformation (cf. section 16. 2\) that preserves the 
field commutation relations. In fact, Z(t) by definition is a squeezed version of A(t) that can be 
obtained from A(t) by passing the latter through a squeezer with the appropriate parameters 
(cf. section IS.2.3j) : in this case the squeezer would have the parameters s = — arccosh(d/^/7) 
and 6 = arg/3. Z(t) satisfies [dZ(t), dZ{t)*\ = dt and the Ito rules for a squeezed field that can 
be generated from the vacuum (the theoretical basis for these manipulations are discussed in 
Appendix [8|) are 



Q T dt. 



Hi n t can be implemented in one arm of a ring cavity with a fully reflecting mirror M and a 
partially transmitting mirror M' with coupling coefficient 7, with Z(t) incident on M'. After 
the interaction, an output field Z ou t(t) is reflected by M' given by 



Zout{t) = u(tyz(t)u(t) 



a 



U(t)*A(t)U{t) 



A 

V7 



U(t)*A{t)*U(t). 



However, the actual output that is of interest is the output Y(t) = U (t)* A(t)U (t) when the 
oscillator interacts directly with the field A(t). To recover Y(t) from Z ou t(t), notice that since 
Q is a quasi-unitary transformation, it has an inverse Q^ 1 which is again quasiunitary. Hence 
Y(t) can be recovered from Z ou t(t) by exploiting the following relation that follows directly 
from the fact that (Zi(t), Z^tff = Q(A l (t),A 1 (t)*) T : 



Y(t) 

Y(ty 



Zout(ty 



That is, Y{t) is the output of a squeezer that implements the quasi-unitary transformation Q~ l 
with Z ou t(t) as its input field. The complete implementation of this linear coupling is shown 
in Figure [161 
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Figure 16: Scheme for implementation of a coupling L = aa+/3a* with a > and a > \/3\. Here 



s = — arccosh(5/y / 7), 6 = arg(/3) and the mirror M' has coupling coefficient 7 = y a 2 — \(3\ 2 . 

6.4 Engineering the interactions between one-dimensional open quantum 
harmonic oscillators 

The second necessary ingredient to synthesizing a general generalized open oscillator according 
to Theorem [2] is to be able to implement a direct interaction Hamiltonian H d given by (|15p 
between one-dimensional harmonic oscillators. The only exception to this, where field-mediated 
interactions suffice, is in the fortuitous instance where Rj k and Lj and Sj, j,k = 1, . . . , n, are 
such that H d = 0. The Hamiltonian H d is essentially the sum of direct interaction Hamiltonians 
between pairs of one-dimensional harmonic oscillators of the form H k i = x^Ckixi (k ^ I) with 
Cki a real matrix. Under the assumption that the time it takes for the light in a ring cavity 
to make a round trip is much faster than the time scales of all processes taking place in the 
ring cavity (i.e., the cavity length should not be too long), it will be sufficient for us to only 
consider how to implement H k i for any two pairs of one-dimensional harmonic oscillators and 
then implementing all of them simultaneously in a network. To this end, let a,j = (pj + iqj)/2 
and a*j = (pj — iqj)/2 for j = k, I, and rewrite H k i as 

Hki = eia* k ai + e[a k dl + e 2 (4af + e* 2 a k ai 

for some complex numbers e\ and The first part Hh = £\a* k ai + e^a^a^ can be simply 
implemented by a beam splitter with a mixing angle = |ei|, $ = — arg(ei) + f, \l/ = — <&, 
and H = (see section I6.2.2p . On the other hand, the second part H kl = e2ata* + e\a k ai 
can be implemented by having the two modes a k and a\ interact in a suitable x nonlinear 
crystal using a classical pump beam of frequency 2oj r and effective pump intensity — 2ie 2 in 
a two-mode squeezing process as described in section 16.1.31 The overall Hamiltonian H k i can 
be achieved by positioning the arms of the two ring cavities (with canonical operators x k and 
xi) to allow their circulating light beams to "overlap" at two points where a beam splitter and 
a nonlinear crystal are placed to implement Hh and iT?>, respectively. An example of this is 
scheme is depicted in Figure [T71 



7 Illustrative synthesis example 



Consider a two degrees of freedom open oscillator G coupled to a single external bosonic 
noise field A(t) given by G = (^4x4, Kx, x T diag(i?i, R 2 )x), with x = {qi,Pi, q2,P2) T , K = 
[ 3/2 1/2* 1 i],Ri= L 2 5 ° 3 5 1, and R 2 = [J ° . 

Let G\ and G 2 be two independent one degree of freedom open oscillators given by G\ = 
{hy.2,K x x\,\x\R\X\) and G 2 = (hx2, K 2 x 2 ,lx^ R 2 x 2 ), with x\ = (<7i,Pi) T , x 2 = fe,P2) r , 
K\ = [ 3/2 i/2 ], and K 2 = [ 1 i ]. Since the scattering matrix for G is an identity matrix, 
it follows from Theorem [2] that G may be constructed as a reducible network given by G = 
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Figure 17: Example implementation of the total direct interaction H^i 
the modes a& and a; of two ring cavities. 



Hl t + Hi, between 



{{Gi, G2}, Hf 2 , G2<G\} with the direct interaction Hamiltonian Hf 2 between G\ and G2 given 
by (cf. (USD) 



zjd 



2i 

1 rp 

2 X2 



1 -xI{k\k x - K^K*)xi 



-1 
3 



Xi. 



This network is depicted in Figure [5j 

In the following we shall illustrate how to build G\ and G2 and how Hf 2 can be implemented 
to synthesize the overall system G. 



7.1 Synthesis of G 1 and G 2 

Let us now consider the synthesis of G\ = (I2x2,KiXi,\x 1 [RiXi). From the discussion in 



section 



2 

0.5 



0.5 

3 



can be realized as a DPA with parameters A = 5 and e = 1 + i, 

while the coupling operator L\ = K\X\ can be realized by the first scheme proposed in section 
6,31 and shown in Figure [15] by the combination of a two-mode squeezer, a beam splitter, and 
an auxiliary cavity mode. Suppose that the coupling coefficient of the mirror M is 72 = 100; 
then the effective pump intensity of the two-mode squeezer is set to be 10 and the beam 
splitter should have a mixing angle of —10 with all other parameters equal to 0. Overall, 
the open oscillator G\ with Hamiltonian H\ = ^x^RiXi and coupling operator L\ can be 
implemented around a ring cavity structure, as shown in Figure [T8l The open oscillator G2 can 
be implemented in a similar way to G\. The Hamiltonian H2 = \x\R2X2 can be implemented 
in the same way as H\ with the choice A = 2 and e = 0. Since e = 0, this means no optical 
crystal and pump beam are required to implement R2, but it suffices to have a cavity that is 
detuned from tj r , the reference frequency in Remark [H by an amount A = 2. The coupling 
operator L2 = 52 + ip2 = 2a2, where a>2 is the annihilation operator /cavity mode of cavity is 
standard and can be implemented simply with a partially transmitting mirror with coupling 
coefficient k = 4, on which an external vacuum noise field A2{t) interacts with the cavity mode 
ci2 to produce an outgoing field Y2(t). The implementation of G2 is shown in Figure [T9l 
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Figure 18: Realization of G\. 




Figure 19: Realization of G<i- 



7.2 Synthesis of Hf 2 



We now consider the implementation of the direct interaction Hamiltonian Hf 2 given by Hf 2 



-x T 



2" L 2 



x\. To proceed, we first note that Hf 2 may be reexpressed in terms of the cavity 

modes a\ and 02 as Hf 2 = —i{a\a2 — a\a 2 ) + 2i(a\a 2 — 0,1(12). Define Hf 2 1 = —i{a\a2 — a±a 2 ) 
and H( 2 2 = 2i(ala 2 -aia 2 ) so that Hf 2 = Hf 21 + H( 2 2 . The first part H( 2 1 = -i(a\a2- a±a 2 ) 
can be simply implemented as a beam splitter with a rotation/mixing angle = — 1 and all 
other parameters set to (cf. section f6 . 2 . 2 1) . On the other hand, the second part Hf 22 = 
2i(a\a 2 — 0,10,2) can be implemented by having the two modes and ai interact in a suitable 



X 



( 2 ) nonlinear crystal using a classical pump beam of frequency uj p = 2u r and effective intensity 



7.3 Complete realization of G = G 2 }, Hf 2 , G 2 < 

The overall two degrees of freedom open oscillator G can now be realized by (i) positioning the 
arms of the two (ring) cavities of G\ and G2 to allow their internal light beams to "overlap" 
at two points where a beam splitter and a nonlinear crystal are placed to implement Hf 2 1 and 
Hf 2 2 , respectively, and (ii) passing the output Y\{t) of G\ as input to G2. This implementation 
is shown in Figure [20j 



8 Conclusions 

In this paper we have developed a network theory for synthesizing arbitrarily complex linear 
dynamical quantum stochastic systems from one degree of freedom open quantum harmonic 
oscillators in a systematic way. We also propose schemes for building the one degree of freedom 
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Figure 20: Realization of G. The block diagram at the top shows how G is realized by a series 
connection of G\ into G2 and a bilinear direct interaction Hf 2 between the canonical operators 
of G\ and Gi- The bottom figure shows the physical implementation of G based on the block 
diagram. 



oscillators and the required interconnections and interactions among them, in the setting of 
quantum optics. 

Together with advances in experimental physics and the availability of high-quality basic 
quantum devices, it is hoped the results of this work will assist in the construction of high- 
performance coherent linear quantum stochastic controllers and linear photonic circuits in the 
laboratory for applications in quantum control and quantum information science. 

Adiabatic elimination of coupled cavity modes 

In this section, we shall derive formulas for two coupled cavity modes in which one of the 
cavities has very fast dynamics compared to the other and can be adiabatically eliminated, 
leaving only the dynamics of the slow cavity mode. The cavities are each coupled to separate 
bosonic fields and are interacting with one another in a classically pumped nonlinear crystal. 
A mathematically rigorous theory for the type of adiabatic elimination/singular perturbation 
that we are interested in here has recently been developed in |33j. 

The two cavity modes will be denoted by a and b, each defined on two distinct copies of 
the Hilbert space Z 2 (Z + ) of square-integrable sequences (Z+ denotes the set of all nonnegative 
integers). Thus the composite Hilbert space for the two cavity modes is H = / 2 (Z + ) ® l 2 (Z, + ). 
The interaction in a nonlinear crystal is given, in some rotating frame, by an interaction 
Hamiltonian H a b of the form Hob = aa*b + (3a*b* + a*ab* + f3*ab for some complex constants a 
and f3. The mode a is coupled to a bosonic field A±, while b is coupled to the bosonic field A2, 
both fields in the vacuum state. The fields A\ and Ai live, respectively, on boson Fock spaces 
T\ and T2, and we denote T = T\ <g> T% . We take a to be the slow mode to be retained and b 
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to be the fast mode to be eliminated. 

We consider a sequence of generalized open oscillators G k = (I,U k \H$), with L< fe ) = 
(y/jia, k^&b) T and H$ = A x a*a + k 2 A 2 b*b + k(aa*b + (3a*b* + a*ab* + (3*ab) each evolving 
according to the unitary Uk satisfying the left H-P QSDE (as opposed to the right H-P QSDE 
in g])): 

dU k (t) = U k (t)(m^dA 1 (t),dA 2 (t)) T -m T (dA 1 (t),dA 2 (t))^ 

_ Imwdt 

Here we are using the left QSDE following the convention used in [33] (see Remark 2 therein) 
so that the interaction picture dynamics of an operator x is given by xit) = Uk(t)xUk(t)* ■ We 
shall use the results of |33j to show, in a similar treatment to section 3.2 therein, that in the 
limit as k — > oo: 

lim sup \\U k (t)*<t>-U(t)*<f>\\ =0 M^eUa^T (19) 

fc->oo <t<T 

for any fixed time T > 0, where Hq is an appropriate Hilbert subspace of % (to be precisely 
specified in the next paragraph) for a limiting unitary U (t) (again as a left H-P QSDE) satisfying 

dU(t) = U{t) (( %2 * 2 + 72 - l) dA 22 + ^a*dAx(t) - SnadAxit)* (20) 
\\z2A 2 -j 2 J 

- (iA 2 - 7 |)" 1 (aa* + P*a)dA 2 (t) 

+ i («*« + Pa*)dA 2 (ty + (iA 1 - —) a*adt 

%2A 2 - 72 v 2 / 

+ (iA 2 - -j) 1 {aa* + (3*a)(a*a + (3a*)dt 

on Hq ® T . Note that (i20l) is a left H-P QSDE corresponding to the right form in section [2] by 
noting that we may write 

%A X - — ) a* a + (iA 2 - — ) [aa* + /3*a)(a*a + (3a*) 



2 J V 2 

= i \A ia *a - A 2 + (72 )2 ( aQ * + + £a*)J " + t\U\ 

with L\ = ^/tTo and L 2 = i^/ 2 (—iA 2 — ^-)~ 1 (a*a + (3a*). As such, it satisfies the H-P 
Condition 1 of [33]. 

Let <po,4>i,... be the standard orthogonal bases of Z 2 (Z + ), i.e., 4>i is an infinite sequence 
(indexed starting from 0) of complex numbers with all zeros except a 1 in the Ith place. First, 
let us specify that Hq = l 2 (%+) <S> C(f>o; this is the subspace of H where the slow dynamics of 
the system will evolve. Next, we define a dense domain V = span{^ ® (pi; j, Z = 0,1,2,...} of 
"H. The strategy is to show that \33\ Assumptions 2-3] are satisfied, from which the desired 
result will follow from [33] Theorem 11]. 

From the definition of H^ , and U k given above, and we can define the operators 
Y,A,B,Gx,G 2 , and W jt = 1,2) in [331 Assumption 1] as follows: Y = (iA 2 - ^)b*b,A = 
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i(aa*b + a*ab* + f3a*b* + (3*ab),B = (iAi - £)a*a,Gi = ^na*,G 2 = 0,F X = 0,F 2 
-y/T^fe*, Wji = 5j[. Then we can define the operators K^ k \ in this assumption as 



K (k) = k 2 Y + kA + B ^ L (k) = + Q (j = 1, 2). 

Let .Po be the projection operator to T~Lq. Let us now address Assumption 2 of [33]. From 
our definition of Hq, it is clear that we have that (a) PqD C T>. Any element of PqT> is of 
the form / ® 0o for some / G span{0;; Z = 0,1,2,...}; therefore, since Y = (iA 2 — ^-)b*b 
and bcpo = 0, we find that (b) YP^d = V d G P. Define the operator Y on V defined by 
F/ ® fa = and Y/ ® 0* = - ^ ® 0j for Z = 1, 2, . . . (Y can then be defined 

to all of T> by linear extension). From the definition of Y and Y, it is easily inspected that 
(cl) YYf = YYf = P\f V/eP, where P\ = I — Pq (i.e., the projection onto the subspace 
of "H complementary to / Hq). Moreover, because of the simple form of Y, it is also readily 
inspected that (c2) Y has an adjoint Y* with a dense domain that contains T>. Since F\ = 0, 
we have that (dl) F?P = on V, while since F 2 */ ® O = y/j^bf ® O = V/ G / 2 (Z + ), we 
also have (d2) F^Pq = on P. Finally, from the expression for A and the orthogonality of the 
bases 0o, 0i, . . ., a little algebra reveals that (e) PoAP d = V d G P. From (a), (b), (cl-c2), 
(dl-d2), and (e), we have now verified that Assumption 2 of [33] is satisfied. 

Finally, let us check that the limiting operator coefficients K, L\ , L 2 , Mi, M 2 , and Nj/. (i, j = 
1, 2) (as operators on 'Ho) of Assumption 3 of [33^ coincide with the corresponding coefficients of 
(1201) . These operator coefficients are defined as K = P (B - AYA)P , L j = P (Gj - AYFj)P , 
M i = -El=iPoW jr (G* r - F r *YA)P , and ity = Y? r =i P Q W jr {F*YF l + 5 rl )P . From these 
definitions and some straightforward algebra, we find that for all / G span{0;; I = 0, 1, 2, . . .} 

K/<S>0o = UiA 1 -^a*a+ (^iA 2 -^y\aa* +p*a)(a*a + (3a*yj f0<l) O , 



*A 2 - Y) 



£1/ <8> <fo = \fl\a*f 

L 2 f <2> 0o = -iy 7 ^ ( 

Mi/® 0o = -v^a/ig^ 

M 2 /®0 O = ^(iA 2 -^) 



and 



Mi/® 00 = /® 00, A r i 2 /®0o = O, iV 2 i/ ®0o = O, 

M2/ ® 00 = 1 \ - 0A / ® 00- 

-72 + i2l\ 2 

Therefore, we see that U(t) may be written as 

(2 22 
(Nji - Sji)dA fl + ^2 MjdA* + ^2 LjdAj + Kdt 
3,1=1 i=i i=i 



Since we have already verified that (120p is a bona fide right-QSDE equation, it now follows 
that Assumption 3 of [33] is satisfied. Now (|19p follows from [331 Theorem 11], and the proof 
is complete. 
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Moreover, we can observe from the derivation above that the coupling of a to ^(t) after 
adiabatic elimination will not change if a is also coupled to other cavities modes 63 , . . . , b m via 
an interaction Hamiltonian of the form YULjip^ji ^ + a *ji a *^j + ctj2a*b* + a* 2 abj), and each 
additional mode may also linearly coupled to distinct bosonic fields As, ... , A m , respectively, 
as long as these other modes are not interacting with b and with one another (this amounts 
to just introducing additional operators Fj,Gj, j > 3, etc.). Moreover, under these conditions 
one can also adiabatically eliminate any of the additional modes, and the only effect will be 
the presence of additional sum terms in U(t) that do not involve b, Ai(t), and A2(t). 

Squeezed white noise calculus 

The purpose of this appendix is to briefly recall results from the theory squeezed white noise 
calculus [34] that are relevant as a basis for some formal calculations presented in section 16.31 
As the theory is quite involved, it is not our intention here to discuss any aspects of it in detail, 
but instead to point the reader to specific results of [34] . 

Let J r (L 2 (M)) denote the usual (symmetric) boson Fock space over the Hilbert space L 2 (]R) 
of complex-valued square integrable functions on R. Let Qj? be the Fock vacuum vector, and 
let ao(/) and a^g) for /, g 6 L 2 {M) be the vacuum creation and annihilation operators on 
J r (L 2 (]R)), respectively. Let n € M and c £ C satisfy n > 0, n > 0, and n(n + 1) > c. The 
parameters n, c characterize the so-called squeezed white noise states w n)C [341 section 2.1] that 
are postulated to satisfy the properties (2.1)-(2.5) therein (see also [211 Chapter 10]). However, 
here we will only be interested in the special case of squeezed states with n, c satisfying the 
constraint n(n + 1) = |c| 2 , as this is the special case of squeezed states that can be generated 
from the vacuum state ojq by an appropriate squeezing Bogoliubov transformation [34} equation 
(2.16)]; see [34 \ Theorem 2.3]. It has been shown that the annihilation and creation operators 
0"n,c{f) and a* nc (g) (f,g G L 2 (]R)) corresponding to such a squeezed states can be concretely 
realized as operators on J 7 (L 2 (1R)) [341 Theorem 2.11 part (b)] and are given in terms of the 
vacuum creation and annihilation operators ao(/) and ao(g), as (this follows from |34| Theorem 
2.3 and equation (3.12)]) 

a n ,c{f) = cosh(s)a (/) + e i6l sinh(s)ao(J/), 
< c (/) = cosh( S )oS(/) + e-^srnh( S )a (J/), 

where J : / 1— >■ /*, s = arctanh( 2 ^ 1 ), and 9 = arg(c). Conversely, we have n = ^ cosh(2s) — ^ 
and c = \e %e sinh(2s). The squeezed white noise state 0Jn,c acts on an operator A affiliated to the 
von Neumann algebra n njC (W(L 2 (M))) // of operators on ^"(L 2 ^)) (here n njC (W(L 2 (]R))) de- 
notes the Gelfand-Naimark-Segal representation of the Weyl C*-algebra W(L 2 (]R)) on J 7 (L 2 (R)) 
corresponding to the state io n ,c, and " denotes the double commutant) as 

cj UjC (A) = (top, Ail?), 

where (•, •) is the complex inner-product on J r (L 2 (R)) (antilinear in the first slot and linear in 
the second). 

Let A(t) = a(x[o,t]) be a vacuum bosonic field, where X[o,t] denotes the indicator function 
for the interval [0,t], and define the squeezed bosonic field A nfi (t) = an,c(X[o,t]) with o njC as 
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defined above. Then A UjC and its adjoint A* n c are related to A and A* by 

A n>c {t) = cosh(s)A{t) + e i9 smh(s)A*(t), (21) 
A* ritC {t) = cosh(s)A*{t) +e~ ie sinh(s) A{t). 

Now, consider an open oscillator whose dynamics are given by the H-P QSDE: 

dU(t) = (-iH + dA(t)*L - L*dA(t) - ^Ldt) U(t), (22) 

where H is the quadratic Hamiltonian of the oscillator and L is the linear coupling operator 
to A(t). By using (|2ip and substituting this into the above QSDE, we may rewrite it in terms 
of the A n ^ c and A* n c as follows: 

dU{t) = (- iH + cL4 n ,c(*) f M - MUA n)C (t) (23) 
- ^(nMM* + (n + l)M* M - c*M 2 - cM*M)d?j U(t), 

where M is a new linear coupling operator given by 

M = cosh(s)L + e id sinh(s)L*. 

As shown in [34] , (|23|) can be interpreted on its own as the unitary evolution of a harmonic 
oscillator and a squeezed bosonic field linearly coupled via the coupling operator M, and this 
defines a quantum Markov process on the oscillator algebra (by projecting to the oscillator 
algebra; see [34^ section 3]). In this interpretation of (2), the squeezed bosonic fields A n ^ c and 
A* n c satisfy the squeezed Ito multiplication rules given by 

dA 2 nc = cdt, dA nfi dA* nc = (n + l)dt, dA* nc dA nfi = n, {dA* nc ) 2 = c*dt, 
dA n)C dt = 0, dA* nc dt = 

that forms a basis for a quantum stochastic calculus for squeezed bosonic fields. A formal 
interpretation of this is that (j23[) defines the evolution of a system coupled to A n ^ c via the 
formal interaction Hamiltonian (see [341 section 3.6]): 

H Int (t) = i(M V * n>c (t) - M* Vn>c (t)), (24) 

where ry njC is a squeezed quantum white noise that can be formally written as r\ n)C = a n)C (5(t))Q 
The connection with the discussion in section [6731 is made by identifying the field Z{t) introduced 
therein with A ntC (t), and rj'(t) with rj n c (t). 

x As is often the case, there is technical caveat in that for mathematical convenience the results of [34] are 
derived on the assumption that H and M are bounded operators on the oscillator Hilbert space. Here we do not 
concern ourselves too much with such detail and assume the optimistic view that these results can be extended 
to unbounded coupling operators M, which are linear combinations of the canonical operators of the harmonic 
oscillator, in view of the fact that the left form (cf. Appendix A) of (|22|l . from which the left form of (|23|) 
can be recovered, still makes sense for a quadratic H and the unbounded operator L associated with M (i.e., 
L — cosh(s)Af — e~ l6 sinh(s)M*) [27]. Moreover, singular interaction Hamiltonians of the form ([24} between 
the unbounded canonical operators of a harmonic oscillator and a vacuum or squeezed quantum white noise 
are physically well motivated and widely used in the physics community. See, e.g., |21l Chapters 5 and 10] and 
related references from [34] section 3.6]. 
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